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Abstract. We consider a broad class of stochastic lattice predator-prey models, whose main 
features are overviewed. In particular, this article aims at drawing a picture of the influence 
of spatial fluctuations, which are not accounted for by the deterministic rate equations, on the 
properties of the stochastic models. Here, we outline the robust scenario obeyed by most of 
the lattice predator-prey models with an interaction a la Lotka-Volterra. We also show how a 
drastically different behavior can emerge as the result of a subtle interplay between long-range 
interactions and a nearest-neighbor exchange process. 

1. Introduction. Understanding biodiversity is a central challenge in modern evolution- 
ary and theoretical biology [May] . In this realm, studying population dynamic models to 
understand and identify the mechanisms allowing for coevolution of competing species 
is an important topic. The latter is classically addressed by considering deterministic 
nonlinear differential equations. Within this approach, the set of equations devised many 
decades ago by Lotka and Volterra [Lot] is certainly a paradigm. These authors consid- 
ered two coupled nonlinear differential equations mimicking the evolution of a two-species 
competing system. Hence, within their model, Lotka and Volterra demonstrated that the 
coexistence of the species always occurs and that the densities of the populations regularly 
oscillate in time. In spite of their popularity, the Lotka-Volterra equations have often been 
criticized as being biologically unrealistic and mathematically unstable [May]. Actually, 
to gain some more realistic and fundamental understanding on population coevolution 
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and biodiversity, it is important to take into account discrete spatial degrees of freedom by 
going beyond the deterministic picture. In this context, various stochastic predator-prey 
models have recently been investigated. In this article, based on Refs. [Mobl,Mob2] where 
more details can be found, we aim at drawing a brief overview of the robust properties 
of these lattice systems with Lotka-Volterra interactions. Also, by considering a simple 
ecological model, we will show how the subtle interplay between the degree of mixing of 
a stochastic system and long-range interaction can give rise to surprising features. 

2. Deterministic rate equations. It is useful to first review some results of the deter- 
ministic approach for the two-species system, where predators (A) and prey (B) interact 
according to the following reactions: A ^ (Z) (death rate /x), A + B A + A (preda- 
tion rate A) and A ^ A + A (reproduction rate a) . Neglecting any spatial variations of 
the densities a{t) and b{t) of predators and prey, respectively, one obtains the classical 
Lotka-Volterra equations: 

d ^ Xa{t)b{t) - ^a{t); b = ab{t) - Xa(t)b{t) . (1) 

A linear stability analysis of these equations show that the densities oscillate around the 
center (neutrally stable) fixed point (ac,6c) = ("'/•^iM/'^) with a characteristic frequency 
/ = y/JIa/2TT. In addition, the existence of a conserved first integral K{t) = \[a{t) + 
b(t)] — cr In a(t) — /i In b(t) implies oscillatory kinetics and coexistence in the whole phase 
portrait, which is characterized by cyclic trajectories. 

As the above deterministic cycles are unstable against any perturbations and the so- 
lutions of Eqs. display amplitudes depending on the initial conditions, which is clearly 
an unrealistic feature, the rate equations are often rendered more realistic by including 
growth-limiting terms. For instance, by assuming that the prey carrying capacity is p, 
one is led to the following equation for species B: 

b = crb{t)[l-b{t)/p]-Xa{t)b{t). (2) 

In this case, the new rate equations exhibit an extinction threshold and have two sta- 
ble fixed points: (i) (0,p), corresponding to a system full of prey and extinction of the 
predators, which is a stable node when A < Ac = fi/p (extinction threshold); and (ii) 
([1 — /i/Ap]cr/A, /i/A), which is associated with the coexistence phase and is either a sta- 
ble node (near the extinction threshold) or a focus (deep in the coexistence phase). The 
existence of a Lyapunov function ensures that these fixed points are actually globally 
stable [May,Mobl,Mob2]. 

3. The stochastic lattice Lotka-Volterra model (and its variants). To probe the 
deterministic approach and gain some understanding on the role of spatial fluctuations 
and correlations, we have studied a stochastic lattice Lotka-Volterra model (SLLVM). The 
latter is formulated in the natural language of reaction-diffusion systems and is defined on 
a d— dimensional lattice, on which the above Lotka-Volterra reactions are implemented. In 
our modeling, we mimic spatial limitation of the resources by assuming that each lattice 
site can be at most occupied by one predator or one prey (site restriction). Monte Carlo 
simulations of the SLLVM show that, in dimensions d > 1, the phase portrait displays 
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Figure 1: (color online) Reproduced from [Mob2]. Snapshots of the time evolution (time 
increases from left to right) of the SLLVM on a 512 x 512 lattice, with rates a = 4.0, /i = 
0.1, A = 2.2. The red, blue and dark dots respectively represent the prey, predators and 
empty sites. Initially the system is homogeneous with densities a(0) = 5(0) = 1/3. 



qualitatively the behavior predicted by the rate equations with growth-limiting term (see 
Fig.l of [Mob2]). In fact, it is found that above some critical threshold both predators 
and prey coexist. The related fixed point is either a node or, for high predation rate, a 
focus which is associated with spiraling flows. In the latter case, as shown in Fig. 
rich spatiotemporal patterns of persistent predator-prey "pursuit and evasion" waves (see 
e.g. Murray's book [May]) develop and translate into erratic damped population density 
oscillations [Mob2]. Obviously these features are not captured by the rate equations 
()1I2() . In finite systems, the quasiperiodic fluctuations appear on a global scale, with 
amplitudes vanishing with the system size. We have analyzed the spatial structure and 
the time evolution of the above complex patterns by computing the stationary two-point 
correlation functions and the power spectrum [Mob2] . This has led us to an understanding 
of the typical size of the clusters reported on Fig. rightmost) and has shown that the 
characteristic frequency of the SLLVM is markedly smaller than the ones predicted by the 
deterministic equations [Mob2] . A completely different picture emerges when the reactive 
fixed point is a node, just above the extinction threshold (Ac) for the predators. Here, one 
observes clouds of predators effectively diffusing in a sea of prey (see Fig. 3 of [Mob2]). At 
the critical value, the system undergoes an active-to-absorbing phase transition. Various 
critical exponents were computed (both in 2D and 3D) and found to be in agreement 
with those of the directed percolation (DP) [Jans]. By deriving an equivalent field-theory 
action from the master equation of the SLLVM and by mapping it (near Ac) onto Reggeon 
field-theory, we provided analytic arguments supporting the numerical indications that 
such a transition is in the DP universality class [Mob2] . Similar DP exponents were found 
numerically in many other stochastic predator-prey models (see e.g. [Ant]). As a further 
probe of the robustness of the above scenario, we have considered different variants of the 
stochastic model, namely the SLVVM supplemented with (i) diffusing predators and prey 
(with same diffusivities) ; (ii) biased diffusion of both species; (iii) predators "following" 
the prey, with biased hoping rates towards neighboring sites occupied by prey. In all these 
cases, the above features were reproduced. Let us also note that for SLLVM without site 
restriction it was shown that predators and prey always coexist in ID and 2D [Wash]. 
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Figure 2: (color online) Reproduced from [Mobl]. Effect of stirring on the phase portrait 
of the 2D stochastic lattice predator-prey with NNN interaction (rates rj = fi = a = 1, 
5 — 10, system size: 256 x 256). ^From left to right: 'stirring' rate V — 0,2, 5. 

4. Mean-field behavior through short-range stirring. We have outlined the robust 
properties of most lattice predator-prey models. Among these characteristics, we have 
seen that deep in the coexistence phase the systems display complex clusters of activities 
at the interfaces of which the dynamics takes place. These features are not affected by the 
diffusion of either predators or prey. On the other hand, it is important to identify if there 
are ecological mechanisms able to render stochastic predator-prey models more tractable 
by means of deterministic equations. This requires to find processes allowing to efficiently 
stir the system by bringing the reactants within the interfaces of the clusters. Arguably, 
the simplest candidate is the process allowing any neighboring sites to exchange their 
content (a prey can avoid an incoming predator): X + Y Y + X, with X,Y E {A, B, 0} 
and X 7^ F. To test the efficiency of this ingredient, we have focused on a model whose 
mean- field (MF) behavior is markedly different from that of Eqs. (|1I2|) . Actually, it is 
natural to split the Lotka-Volterra predation reaction by introducing two independent 
time scales and the new processes: (i) reproduction of the predators in the vicinity of 
a prey {favorable environment), according to the next-nearest-neighbor (NNN) reaction 
A-\-(Z)-\- B A-\-A-\-B (with rate proportional to 6); and (ii) consumption of a prey by a 
neighboring predator, as A + B ^ (Z) + B (with rate proportional to r/). The MF equations 
for this model, which can be viewed as comprising also a nearest-neighbor reproduction 
reaction {A (Z) ^ A -\- A) taking place on a much longer time scale, read [Mobl] 

d = Sa{t)b{t)[l - a{t) - b{t)] ~ iia{t); b ^ ab{t)[l - a{t) - b{t)] - r]a{t)b{t). (3) 

In stark contrast to Eqs. pi2(l . these equations admit two stable fixed points above a 
given threshold 5c'. one is (always) a node and corresponds to an absorbing steady state 
(system full of prey), while the other fixed point (either a node or a focus) is reactive 
(coexistence of prey, with a density < 1/2, and predators) . Hence, according to MF theory, 
for 5 > 5c this model undergoes a first- order phase transition. To check when these MF 
predictions actually apply to the stochastic lattice model, the latter has been numerically 
simulated in the presence of the exchange process (with stirring rate proportional to 
V) [Mobl]. It turns out that there is a subtle interplay between the NNN interaction 
and the exchange process. As illustrated in Fig. [3 the latter can be summarized as 
follows: (a) For vanishing mixing {T> small compared to the other rates), fluctuations 
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have a drastic effect and invalidate the MF picture in dimensions 1 < d < 4, where 
the system undergoes an active-to-absorbing state transition belonging again to the DP 
universality class, (b) When one allows for short-range particle exchange {V > 0), the 
phase portrait flows change dramatically (Fig. 12 center), (c) When the exchange process 
becomes sufficiently fast (e.g. V ^ S) a, reactive fixed point is also available (in any 
dimension), as demonstrated in Fig. El (right). In this case, the system undergoes a first- 
order phase transition as predicted by the MF theory and, for 'fast' stirring, the latter 
becomes very accurate. We have also found that the stable reactive fixed point is either a 
node or a focus. In the latter case, the coexistence phase is again characterized by moving 
activity fronts but, as the system is more mixed, these clusters appear less prominent than 
in Fig.[T][Mobl]. 

5. Conclusion. We have outlined the robust features of a class of stochastic lattice 
models with Lotka-Volterra interactions and discussed how their deterministic (mean- 
field) descriptions are altered by spatial fluctuations and correlations. Wc have also shown 
that the rate equations can aptly describe the dynamics of a stochastic model in the 
presence of an efficient short-range exchange process. We have illustrated this point by 
considering a system with NNN interaction which exhibits either a first- or second-order 
phase transition, depending on the stirring rate. 
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